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Abstract 

Spin models of adaptive networks have been shown to exhibit self-organized critical behavior from simple, 
biologically inspired and locally defined adaptation rules, and have been used to argue for similar adaptive 
mechanisms in the brain [1]. However, the choice of spin- type node states limits the biological plausibility 
of these models. Subsequent experimental studies on real cortical cultures have revealed avalanche-like 
propagation of activity with statistical properties that support the hypothesis of neural dynamics near 
criticality [5] . We here translate the spin model to a network of biologically plausible threshold nodes with 
Boolean states and adapt the original correlation-based rewiring algorithm. As a result the new model 
allows to quantify neuronal avalanches in the evolved networks, where power-law scaling exponents of 
avalanche sizes and durations are compatible both with relations predicted from universal scaling theory 
as well as results from observed activity avalanches in real cortical cultures. 

Introduction 

Information processing by a network of dynamical elements is a delicate matter: Avalanches of activity 
can die out if the network is not connected enough or if the elements are not sensitive enough; on the 
other hand, activity avalanches can grow and spread over the entire network and override information 
processing as observed in epilepsy. Therefore, it has long been argued that neural networks have to 
establish and maintain a certain intermediate level of activity in order to keep away from the regimes of 
chaos and silence [TJ|3H5]. Similar ideas were also formulated in the context of genetic networks where 
Kauffman postulated that information processing in these evolved biochemical networks would be optimal 
near the "edge of chaos" , or criticality of the dynamical percolation transition of such networks [5] . 

In the wake of self-organized criticality (SOC), it was asked if also neural systems were self-organized 
to some form of criticality [7J. An early example of a SOC model that had been adapted to be appli- 
cable to neural networks is the model by Eurich et al. [8] . They show that their model of the random 
neighbor Olami-Fcdcr-Christcnscn universality class exhibits (subject to one critical coupling parameter) 
distributions of avalanche sizes and durations which they postulate could also occur in neural systems. 

Another early example of a model for self-organized critical neural networks [TJG3 drew on an alter- 
native approach to self-organized criticality based in dynamical networks |10j . Here networks are able 
to self-regulate towards and maintain a critical system state, via simple local rewiring rules which are 
plausible in the biological context. 

After these first model approaches, indeed strong evidence for criticality in neural systems has been 
found in terms of spatio-temporal activity avalanches, first in the seminal work of Bcggs and Plenz [2J. 
Much further experimental evidence has been found since, which we will briefly review below. These 
experimental findings sparked intense research on dynamical models for criticality and avalanche dynamics 
in neural networks, which we also give a brief overview below. While most models emphasized biological 
and neurophysiological detail, our path here is different: The purpose of our paper is to pick up the 
thread of the early self-organized critical neural network model [T| and test its applicability in the light of 
experimental data. We would like to keep the simplicity of the first spin model in the light of statistical 
physics, while lifting the drawback of a spin formulation w.r.t. the biological system We will study 
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an improved model and show that it adapts to criticality exhibiting avalanche statistics that compare 
well with experimental data without the need for parameter tuning. 

Avalanche dynamics in neuronal systems 

Let us first briefly review the experimental studies on neuronal avalanche dynamics. In 2003, Beggs and 
Plenz published their findings about a novel mode of activity in neocortical neuron circuits [2J. During 
in-vitro experiments with cortex slice cultures of the rat, they found evidence of spontaneous bursts 
and avalanche-like propagation of activity followed by silent periods of various lengths. The observed 
power-law distribution of event sizes indicates that the neuronal network is maintained in a critical state. 
Also, the spatio-temporal patterns of the avalanches are stable and precise over many hours and robust 
against external perturbations |12| . which indicates they might play a central role for brain functions 
like information storage and processing. Neuronal avalanches have also been found during developmental 
stages of in-vitro cortex slice cultures from newborn rats [13j , as well as in cultures of dissociated neurons 
in different kinds of networks, like rat hippocampal neurons and leech ganglia |14j or rat embryos |15j . 

Aside from these in-vitro experiments, extensive studies in-vivo have since been conducted. The 
emergence of spontaneous neuronal avalanches has been shown in anaesthesized rats during cortical 
development |16) as well as in awake rhesus monkeys during ongoing cortical synchronization |17j . 

The biological relevance of the avalanche-like propagation of activity in conjunction with a critical 
state of the neuronal network has been emphasized in several works recently. Such network activity 
has proven to be optimal for maximum dynamical range [18 ] ll9j. maximal information capacity and 
transmission capability [20] as well as for a maximal variability of phase synchronization |21j . Most 
recently, experimental evidence for universality of critical dynamics has been found in neuronal avalanche 
data [22] and formally linked to universal scaling theory |23j . 

These experimental studies with their rich phenomenology sparked a large number of theoretical 
studies and models for criticality and self-organization in neural networks, ranging from simple toy models 
to detailed representations of biological functions. 

A variety of models have been constructed close to biological details at the neuron level, and e.g. 
use threshold firing dynamics and activity-dependent plasticity of synaptic couplings as the basis for 
self-organization. While some models feature synaptic facilitation following a firing event [24H26] , others 
use synaptic depression as the main driving force towards criticality |27|I28|. It has been shown that 
anti-Hcbbian evolution is generally capable of creating a dynamically critical network when the anti- 
Hcbbian rule affects only symmetric components of the connectivity matrix, while the anti-symmetric 
component remains as an independent degree of freedom utilizable for e.g. learning tasks |29| . Also 
the biological plausibility of activity-dependent synaptic plasticity for adaptive self-organized critical 
networks has been stressed [30]. Recently, correlations of subsequent firing events again came into focus 
as a synaptic facilitation criterion [31] . The biological relevance of the critical state in neural networks for 
a brain function as learning has further been underlined [32] . Most recently, the temporal organization 
of neuronal avalanches in real cortical networks has been linked to the existence of alternating states 
of high vs. low activity in the network as well as to a balance of excitation and inhibition in a critical 
network [33] . 

While the proposed organization mechanisms strongly differ between the individual models, the re- 
sulting evolved networks tend to be part of the same fundamental universality classes, exhibiting e.g. 
avalanche statistics in a similar way as in the experimental data, as power-law distributions at exponents 
of —3/2 for avalanche sizes. With the recent, more detailed models in mind, we are especially interested 
in the underlying universality of self-organization. 

Considering the enormous interest in neuronal self-organization, we here come back to our older spin 
models (HQS] and develop a new basic mechanism in the light of improved biological plausibility of these 
models. Two main aspects have to be addressed: In a previous work [llj . we have pointed out the 
problems of a spin-type description in the context of biological, esp. neural, networks, and discussed a 
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new basic model featuring Boolean state nodes and a redefined activation threshold function. We will use 
the Boolean state model in the following. As a major advantage, activity avalanches intrinsically occur in 
this type of network, whereas spin networks typically exhibit continuous fluctuations with no avalanches 
directly visible. However, the correlation-based rewiring mechanism of e.g. pQ will no longer work when 
inactive nodes are now represented by '0' instead of '— 1'. A solution will be presented below. The second 
aspect to be reviewed is topological invariance of the algorithm. The older, correlation-based mechanism 
of self-organization [T] relies on randomly selecting neighboring sites on a lattice for the rewiring processes. 
On a lattice, the number of possible neighbors is strictly limited, but on a large random network near 
critical connectivity, there are far more unconnected pairs of nodes than there are connected pairs. Thus, 
randomly selecting pairs of nodes for rewiring would introduce a strong bias towards connecting nodes 
which were previously unconnected. This results in a strong increase of connectivity and makes a self- 
organized adaptation towards a critical state practically impossible. If we want to overcome the restriction 
of e.g. a confined lattice topology in order to improve biological applicability of the model, we have to 
adapt the rewiring mechanism such that this bias no longer exists. 

Of course, the new model shall inherit several important features from the older spin models, which al- 
ready underline the applicability to biological networks: in particular, it must be capable of self-regulation 
towards a critical state despite being simplified to the most minimal model possible. The organization 
process however should be based on a simple, yet biologically plausible rewiring rule, which only uses 
local information accessible to individual nodes like pre- and post-synaptic activity and correlation of 
such activity. 



Materials and Methods 
Adaptive network evolution 

Consider a randomly connected threshold network of the type discussed in [IT] . The network consists of 
N nodes of Boolean states <7j € {0, 1} which can be linked by asymmetric directed couplings cy = ±1. 
Node pairs which are not linked have their coupling set to = 0; and links may exist between any two 
nodes, so there is no underlying spatial topology in this model. Let K denote the average connectivity 
of the network, i.e. the number of in- links averaged over all N nodes. 

All nodes are updated synchronously in discrete time steps via a simple threshold function of their 
input signals with a little thermal noise introduced by the inverse temperature /3, in the same way as 
in the model of Bornholdt and Roehl [TJ, but now with an input shift of —0.5 in the Glauber update, 
representing the new Oo function defined in 

Probst + 1) = 1] = gp(Mt)) 

Prob[ CTi (t + l) = 0] = 1 (1) 

with 

N 

/<(*) = I>i*i(*)-e< (2) 

and 

For the simplicity of our model, we first assume that all nodes have an identical activation threshold of 
&i = 0. 

The adaptation algorithm is now constructed in the following way: We start the network at an 
arbitrary initial connectivity Ki n i and do parallel synchronous updates on all nodes according to eq. 
([T]). All activity in this model originates from small perturbations by thermal noise, leading to activity 
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avalanches of various sizes. In this case we set the inverse temperature to (3 = 5. On a larger time 
scale, i.e. after r = 100 updates, a topology rewiring is introduced at one randomly chosen, single 
node. The new element in our approach is to test whether the addition or the removal of one random 
in-link at the selected node will increase the average dynamical correlation to all existing inputs of that 
node. By selecting only one single node for this procedure, we effectively diminish the bias of selecting 
mostly unconnected node pairs - but retain the biologically inspired idea for a Hcbbian, correlation-based 
rewiring mechanism on a local level. 

Now, we have to define what is meant by dynamical correlation in this case. We here use the Pearson 
correlation coefficient to first determine the correlation between a node i and one of its inputs j over the 
preceding r time steps: 



where Si and Sj in the denominator depicts the standard deviation of the state of nodes i and j respec- 
tively. In case one or both of the nodes are frozen in their state (i.e. yield a standard deviation of 0), 
the Pearson correlation coefficient would not be defined, we will assume a correlation of Cij = 0. Note 
that we always use the state of node i at one time step later than node j, thereby taking into account the 
signal transmission time of one time step from one node to the next one. Finally, we define the average 
input correlation C™ 9 of node i as 



where k% is the in-degree of node i. The factor |cy| ensures that correlations are only measured where 
links arc present between the nodes. For nodes without any in-links (ki = 0) we define Cf 1 ' 9 := 0. 
In detail, the adaptive rewiring is now performed in the following steps: 

1. Select a random node i at which an eventual rewiring will take place. 

2. Run network updates for r = 100 simulation time steps and measure the average input correlation 
C™ 9 of node i. 

3. With equal probability, either 

(a) insert an additional in-link of random weight a x = ±1 at node i from a random, previously 
unlinked node x, or 

(b) remove one of the existing in-links at node i. 

4. Again run r = 100 network updates and measure the new Cf 9 of node i after the local topology 



5. If Cf 9 has increased after the insertion or removal of the in-link, the topology change from step 3 
is finally accepted; otherwise, the rewiring is reverted. 

6. Run r = 100 network updates to allow for an eventual transient period to be completed prior to 
the next rewiring process. Iterate from step 1. 

Note that the exact choice of r is not critical, but is chosen as r = 100 here to ensure time scale seperation 
of at least two orders of magnitude between node dynamics (fast) and topology changes (slow). 

It is also worth noting that this model - in the same way as the earlier work pQ - uses locally available 
information at synapse level and takes into account both pre- and post-synaptic activity. This is a 
fundamental difference to approaches discussed e.g. in [24], [25] or [27], where only pre-synaptic activity 
determines changes to the coupling weights. 




(4) 




(5) 



change. 
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In order to obtain an indication of the current dynamical regime of the network (i.e. whether the 
network is sub- or super-critical or close to the critical point), we continuously measure a branching 
parameter based on potential damage spreading in the network. This is realized by counting, for each 
individual node i, the number of descendant nodes which would possibly change their states at time step 
t + 1 if the state of node i was changed at the present time t. Here, both the present states (on or off) 
of node i and its descendants as well as the nature of their respective links (activating or inhibiting) arc 
taken into account. The obtained number of descendant nodes prone to damage spreading (and thus also 
to signal propagation) is then averaged over the entire network, resulting in the branching parameter A. 
This allows to estimate (based on the current network configuration) whether the network is sub-critical 
(A < 1) or super-critical (A > 1). 

For the later analysis of e.g. avalanche statistics in the evolved critical networks, we export snapshots 
of the network structure whenever the branching parameter is close to one (typically when |A — 1| < 0.01). 



Avalanche analysis 

For a detailed analysis of avalanche statistics, we use the snapshots of near-critical (per branching param- 
eter estimation) network structures from the adaptation run explained above. Key observables are the 
avalanche size S, i.e. the total number of nodes which become active at least once during one avalanche, 
and the avalanche duration T, i.e. the number of simulation time steps from the start (first node active) 
to the termination (no more nodes active) of an avalanche. To obtain those, the network is now run in a 
deterministic mode with j3 — > oo. The update function from (fT]) then simplifies to 

*i(t + i) = e (fi(t)) (6) 

with a redefined threshold function (nodes only become active with activating, non-zero input) 

& ^={l: ill p> 

and the usual input function 

N 

f i (t)=Y,c i w(t)-0 i (8) 

where activation thresholds 0i are set to first. With parallel updates, any network activity would 
eventually end up in either a fixed point or limit cycle attractor, but not necessarily at the fixed point 
"all nodes off" , terminating an avalanche. Therefore, we introduce an exhaust time parameter 57 which 
can be biologically interpreted as an effect of depleting neuro-transmitter substances at active synapses. 
In each time step, every node will increase its individual activation threshold 9i to 1 with a probability 
corresponding to its own activity (i.e. number of time steps where <7j = 1) over the last Q time steps. This 
turns out to be sufficient in eventually stepping out of a periodic attractor and terminating the avalanche. 
Whenever one avalanche is terminated (all nodes off), we will start a new one by randomly activating one 
single node and continue with the parallel updates. We constantly keep track of cumulative avalanche size 
and duration distributions, as well the average size (S)(T) of avalanches that have a certain duration T. 
From universal scaling theory |23j we expect the following power-law scaling relations in case of critical 
networks: 

fcml(S) ~ S~ T ™' (9) 

fc m i(T)~T- a ^' (10) 

(S)(T) r*j r p l / auz (11) 

where the exponents fulfill 

£W = J_ ( 12 ) 
T cm i avz 
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Results 

Adaptive network evolution 

In the following, we will have a look at different obscrvablcs during numerical simulations of network 
evolution in the model presented above. Key figures include the average connectivity K and the branching 
parameter A. Both arc closely linked to and influenced by the ratio of activating links p, which is simply 
the fraction of positive couplings Cjj = +1 among all existing (non-zero) links. 

The upper part in Figure Q] shows an exemplary run of the topology evolution algorithm, where we 
start with completely isolated nodes without any links. Trivially, the "network" is subcritical at this 
stage, which can be seen from the branching parameter which is far below 1. As soon as rewiring takes 
place, the network starts to insert new links, obviously because these links enable the nodes to pass signals 
and subsequently act in a correlated way. With increasing connectivity, also the branching parameter 
rises, indicating that perturbations start to spread from their origin to other nodes. When the branching 
parameter approaches 1, indicating that the network reaches a critical state, the insertion of new links 
is cut back. The processes of insertion and depletion of links tend to be balanced against each other, 
regulating the network close to criticality. 

If we, on the other hand, start with a randomly interconnected network at a higher connectivity like 
Ki n i = 4 as shown in the lower part of Figure [1] we find the network in the supercritical regime (A > 1) 
at the beginning. When above the critical threshold, many nodes will show chaotic activity with very 
low average correlation to their respective inputs. The rewiring algorithm reacts in the appropriate way 
by deleting links from the network, until the branching parameter approaches 1. 

In both examples above, the evolution of the ratio of activating links p (which tends towards 1) shows, 
that the rewiring algorithm in general favors the insertion of activating links and vice versa the deletion 
of inhibitory couplings. This appears indeed very plausible if we remind ourselves that the rewiring 
mechanism optimizes the inputs of a node towards high correlation on average. Also, nodes will only 
switch to active state <7j = 1 if they get an overall positive input. As we had replaced spins by Boolean 
state nodes, this can only happen via activating links - and that is why correlations mainly originate from 
positive couplings in our model. As a result, we observe the connectivity evolving towards one in-link 
per node, with the ratio of positive links also trending towards one. 

For a richer pattern complexity, we might later want to introduce a second mechanism which balances 
out positive and negative links better, and with a first approach we can already test how the rewiring 
strategy would fit to that situation: if, after each rewiring step, we change the sign of single random 
links as necessary to obtain a ratio of e.g. 80% activating links (i.e. p = 0.8), keeping the large majority 
of present links unchanged, the branching parameter will again stabilize close to the critical transition, 
while the connectivity is maintained at a higher value. Figure [2] shows that the self-organization behavior 
is again independent from the initial conditions. This result does not depend on the exact choice of the 
activating links ratio p; similar plots can easily be obtained for a large range starting at p = 0.5, where 
the connectivity will subsequently evolve towards a value slightly below K = 2, which is the value we 
would expect as the critical connectivity for a randomly wired network with balanced link ratio according 
to the calculations made for the basic network model 

Avalanche properties 

In addition to the branching parameter measurement, we also take a look at some dynamical properties 
of the evolved networks to further verify their criticality. As stated in the introduction section, we are 
especially interested in eventually resulting activity avalanches on the networks. Figure [3] shows the 
distributions of avalanche characterics like size and duration, as well as indications that the resulting 
avalanches on the evolved networks show universal properties across scales. (A): Avalanche size S shows 
a power-law scaling f cm i(S) ~ S~ Tcml almost up to network size with an exponent T crn i = 0.5 ± 0.1 in the 
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Figure 1. Exemplary run of the network self-organization algorithm. Regardless of initial 
connectivity and dynamical regime, the network evolves to a critical configuration. Top: when starting 
with completely isolated nodes the "network" is obviously subcritical and links will be inserted. Thus, 
both the connectivity (red) and the branching parameter (blue) as an indicator of network criticality 
increase. The network approaches a critical state where the branching parameter stabilizes close to one. 
Bottom: with higher initial connectivity, the network is supercritical at first. Links are removed from 
the network while the branching parameter approaches the critical value of one. As the 
self-organization algorithm is constructed to optimize activity correlations between linked nodes, the 
ratio of activating links (green) slowly increases in both cases. 

cumulative distribution, corresponding to a probability density exponent ofr = 1.5±0.1. (B): Similarly, 
avalanche durations T are power-law distributed as well up to a duration of approximately 70 time steps, 
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Figure 2. Exemplary run with fixed ratio of activating links. If the ratio of activating links 
(green) is kept fixed (here: p = 0.8; i.e. 80% activating links) in order to keep some inhibiting links 
within the network, the connectivity (red) evolves to a higher value. Still, the branching parameter 
(green) is maintained close to the critical value of one. Top: starting with isolated nodes (subcritical). 
Bottom: starting at supercritical connectivity. 



according to f cm i(T) ~ T~° cml with an exponent of a cm i = 0.9 ± 0.1, i.e. a = 1.9 ± 0.1. As discussed 
above, plain power-law distributions can originate from several mechanisms and must not seen as clear 
evidence of criticality. To obtain a third exponent, we have also measured the average avalanche sizes 
(S)(T) as a function of avalanche duration T. It becomes clear from (C) that for durations of approx- 
imately 70 time steps and more, the avalanches begin to span most of the system size, which explains 
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the cutoff position in the avalanche duration scaling (B). Up to that point, we find a power-law scaling 
(S)(T) ~ j 1 !/ ^ with an exponent of = f .8 ± 0.1. These exponents are both in line with experimen- 
tal results [35] and fulfill the exponent relation = = predicted for a critical system |23) . 
(D): Finally, we find that avalanche profiles (i.e. average scaled size as a function of scaled avalanche 
duration) of avalanches of different durations T approximately collapse onto a universal shape. In other 
words, avalanches are self-similar across different scales, meaning e.g. each small avalanche could be part 
of or could possibly trigger larger avalanches. 





1 10 100 1000 0.2 0.4 0.6 0.8 1 

avalanche duration T fraction of avalanche duration (t/T) 



Figure 3. Critical exponents and collapse of avalanche profiles. Results measured from 10 5 
avalanches on 10 different evolved sample networks of TV = 1024 nodes, with an exhaust time parameter 
il = 1000. Similar results were gained with ft = 100 or O = 10000, this choice has no significant effect 
on the scaling exponents. A: Cumulative distribution of avalanche sizes shows a power-law scaling 
fcmi(S) ~ S~ T " nl with exponent r cm / ~ 0.5. B: Cumulative distribution of avalanche sizes shows a 
power-law scaling f cm i{T) ~ T~ a ° ml with exponent a cm i « 0.9. C: Average size (S) of avalanches of 
duration T shows a power-law increase corresponding to (S)(T) <~ r p^l< TVZ with an exponent of 

m 1.8. Note that the exponents r CTO ;, a cm i, -^z fulfill the relation = which is expected for 
a critical system. D: Profiles of avalanches, i.e. average scaled size as a function of scaled avalanche 
duration, of different durations (shown here for T = 13, 26, 39, 52) approximately collapse onto a 
universal shape. 
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Variations in activation thresholds and response to external perturbation 

In the above simulations, the activation thresholds of all nodes were strictly set to @i = for maxi- 
mum model simplicity. However, a neuron might as well need higher input to become active. Figure 2] 
demonstrates that the proposed adaptation algorithm similarly works well on networks of nodes with an 
increased activation threshold of e.g. Oj = 1. According to the update rule ([T]), now at least two positive 
inputs are necessary to activate a single node. As the rewiring algorithm is based on propagation of 
thermal noise signals, the inverse temperature (3 needs to be selected at a lower value than before. (As 
a general rule, (3 should be selected in a range where thermal activation of nodes occurs at a low rate, 
such that avalanches can be triggered, but are not dominated by noise.) The simulation is now started 
at an average connectivity of K = 7 which is still sub-critical in this case (branching parameter low). 
In a similar way as shown above, the network adapts by inserting new links and increasing K, thereby 
also increasing the average branching parameter. Again, this process is stopped after the supercritical 
regime is reached. While the system does not approach to a phase transition as nicely as shown above 
for activation thresholds of zero (in fact the branching fluctuates much more around the target value of 
one), the general tendency remains: the rewiring mechanism reacts properly before the network drifts 
too far off from criticality. The connectivity also fluctuates more, but stabilizes on a level around K = 9. 
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Figure 4. Exemplary run with higher activation thresholds Oj = 1. When activation 
thresholds are increased, a node needs more than one excitatory input to become active itself. Thus, 
higher overall connectivity is needed to allow signal propagation on a critical level. The adaptation 
process responds accordingly and maintains a connectivity around K = 9 while the branching 
parameter shows larger fluctuations between sub- and supercritical states, but in general is kept on a 
moderate level and does not diverge with increasing connectivity. 

In their in-vitro experiments, Beggs and Plenz further demonstrate that cortical networks can self- 
regulate in response to external influences, retaining their functionality while avalanche-like dynamics 
persist - for example after neuronal excitability has been increased by adding stimulant substances to 
the cultures. 

To reproduce such behavior in our model, we can include variations in the activation thresholds Oj 
of the individual nodes. Assume we start our network evolution algorithm with a moderately connected, 
but subcritical network, where all nodes have an activation threshold of Oj = 1. 
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Figure 5. Rewiring response to a sudden decrease of activation thresholds. If we first set all 
node activation thresholds to 0^ = 1, the connectivity (red) must first evolve to a higher value to allow 
propagation of activity within the network. When all activation thresholds are suddenly reduced (to 
mimic an external influence of neuronal excitability by stimulant substances) at the same time step, the 
network properly responds to the new situation and reduces connectivity to decrease excitability back 
to a critical level. 

Figure [5] shows that the network first behaves in the same way as demonstrated in Figure [4] for 
activation thresholds 0j = 1. At one time step in the center of Figure [5l we suddenly reset all nodes 
to an activation threshold of 0.; = 0, simulating the addition of a stimulant. As we can expect, this 
immediately puts the network into a supercritical, chaotic state. This is reflected by the branching 
parameter, which now constantly stays above one and does not fluctuate below anymore. It is clearly 
visible that the rewiring mechanism promptly reacts and drives back the connectivity, until the branching 
parameter slowly converges towards one again. A similar behavior is also found if thresholds 0^ are not 
changed all at once, but gradually during network evolution. 

Discussion 

We have seen that already very minimalistic binary neural network models arc capable of self-organized 
critical behavior. While older models show some drawbacks regarding biological plausibility originating 
directly from their nature as spin networks, we here presented a possible transition to a self-organized 
critical, randomly wired network of Boolean node states with emerging dynamical patterns, namely 
activity avalanches, reminiscent of activity modes as observed in real neuronal systems. This is possible 
via a simple, locally realizable, rewiring mechanism which uses activity correlation as its regulation 
criterion, thus retaining the biologically inspired rewiring basis from the spin version. While it is obvious 
that there are far more details involved in self-organization of real neuronal networks - some of which 
are reflected in other existing models - we here observe a fundamental organization mechanism leading 
to a critical system that may serve as the simplest representative of a "neural SOC universality class" . It 
exhibits statistical properties like avalanche size and duration distributions with scaling exponents that 
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compare well neuronal avalanches found in biological experiments. 
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